Optimal trajectories for efficient atomic transport without final excitation 
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We design optimal harmonic-trap trajectories to transport cold atoms without final excitation, 
combining an inverse engineering techniqe based on Lewis- Riesenfeld invariants with optimal control 
theory. Since actual traps are not really harmonic, we keep the relative displacement between the 
center of mass and the trap center bounded. Under this constraint, optimal protocols are found 
according to different physical criteria. The minimum time solution has a "bang-bang" form, and the 
minimum displacement solution is of "bang-off-bang" form. The optimal trajectories for minimizing 
the transient energy are also discussed. 
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I. INTRODUCTION 

Efficient transport of ultracold atoms and ions by mov- 
ing the confining trap is an important goal in atomic 
physics with applications to basic science, metrol- 

ogy, and quantum information processing. A sufficiently 
slow, adiabatic motion is a simple way to transport the 
atoms without excitations or losses 0-|al- However, the 
long time required may become impractical, e^. if a 
fast quantum-information operation is required P, Q ; or 
counterproductive, because of the accumulation of per- 
turbations. Motivated by this, several theoretical and 
experimental investigations have been devoted to mak- 
ing atomic transport fast, and simultaneously faithful to 
the ideal result of adiabatic transport [9l-[l3| . 

These works on transport share concepts and tech- 
niques with other operations in which a "shortcut to adi- 
abaticity" is to be found, e.g. in expansion or compres- 
sions [nl, [T5I - I23I I , rotations [23| , and internal state pop- 
ulation transfer [25l - l3l| . Several approaches have been 
proposed, including counter-diabatic [25; 27 [or, equiva- 



30|, optimal 
and in- 
Lewis-Riesenfeld invariants 



lently, transitionless driving algorithms 
control theory |19|], "fast- forward" scaling 
verse engineering based on 

In essence, the invariant-based inverse engineering 
method relies on designing the Hamiltonian evolution so 
that the eigenvectors of corresponding invariants of mo- 
tion become at initial and final times equal to the instan- 
taneous eigenvectors of the Hamiltonian. This method 
provides in fact families of paths ^S^] which satisfy the 
initial and final boundary conditions, and thus guarantee 
the fast transitionless evolution, ideally in an arbitrarily 
short time. Given this freedom, it is natural to combine 
the invariant-based inverse method and optimal control 
theory to optimize the trajectory according to different 
physical criteria or operational constraints. For example, 
the time-dependent frequency of a harmonic trap expan- 
sion can be optimized with respect to time or to transient 



excitation energy, with a restriction of the allowed tran- 
sient frequencies (33l - [35l |. 

In this paper, we apply the invariant-based method 
complemented by optimal control theory to find opti- 
mal trajectories for fast atomic transport on harmonic 
traps without final vibrational excitation. Since actual 
traps are not really harmonic, we keep, as an imposed 
constraint, the relative displacement between the center 
of mass and the trap center bounded. We then opti- 
mize the trajectories according to different physical cri- 
teria: time minimization, (time averaged) displacement 
minimization, and (time averaged) transient energy min- 
imization. 



II. INVARIANT-BASED INVERSE 
ENGINEERING METHOD 

We consider here the harmonic transport described by 
the time-dependent Hamiltonian 



H{t) 



P 

2m 



1 



(1) 



where q and p are the position and momentum operators, 
cjQ is the constant harmonic frequency of the potential, 
and qo(t) is the position of the center of the harmonic 
trap. The corresponding quadratic-in-momcntum Lewis- 
Riesenfeld invariant [32] has the form [36, 37] (up to an 
arbitrary multiplicative constant) 

Ht) ^^{P- mq,f + ]^'mujl[q - q,{t)Y , (2) 

where the functions qc{t) must satisfy the auxiliary equa- 
tion 



9c + i^liqc - qo) = 0, 

to guarantee the invariant condition 
dl(t) _ dl{t) 1 



dt 



dt 



I{t),H{t)]=0. 



(3) 



(4) 



2 



Eq. ^ is simply Newton's equation for a classical par- 
ticle in the moving harmonic potential. 

An arbitrary solution of the time-dependent 
Schrodinger equation ihdt^{q,t) — H{t)'i'{q,t), may be 
written in terms of "transport modes" e*"" V'n('Z, i), 



(5) 



where n — 0,1,..., c„ are time-independent coefhcients, 
ipniQ,t) are the orthonormal eigenvectors of the invari- 
ant I{t) satisfying I{t)ipniq,t) = '^ni^niqjt), with real 
time-independent A„, and the Lewis- Riesenfeld phase is 
defined as 

""^^^ " ^ i ("^"^^'^1'^^ ~ (6) 
For the harmonic trap considered here (37l |. 



1 



(2"7i!)i/2 \ nh J 



X exp 



I — -. — ]H„ 



exp 



mujo 



mujQ 
2h 



1/2 



{q - qc) 



(7) 



i.e., qc is the center of mass of the transport modes. Sub- 
stituting Eq. dH) into Eq. ©, 



1 



dt' X 



mq^ 



(8) 



where A„ = En — {n + l/2)hujQ. The instantaneous av- 
erage energy for a transport mode can be obtained from 
Eqs. ((H) and ©, 

(V'nWI-ffWI^nW) = [n + 1/2) + E, + Ep, (9) 

where the first, "internal" contribution remains constant 
for each n, Ec = mq^/2, and Ep = ^mujQ{qc — qo)^ has 
the form of a potential energy for a classical particle. The 
instantaneous average potential energy can be written as 



{Vit)) = ^{n + l/2) + Ep. 



(10) 



Suppose that the harmonic trap is displaced from 
<Zo(0) = to qo{tf) = d in a time tf. The trajectory 
qo{t) of the trap can be inverse engineered by designing 
first an appropriate classical trajectory qc{t). To avoid 
vibrational excitation at the final time we impose the 
conditions 



(7,(0) = 0; qM = 0; <j-,(0) = 0, 
qc{tf) = d; qc{tf) = 0; qdtf) = 0, 

which, along with Eq. ([3]) , imply also 
<Zo(0) =0; qo{tf)^d. 



(11) 
(12) 



(13) 



The above boundary conditions guarantee the commu- 
tativity of I{t) and H{t) at t = and t — tf, that is, 
the transport modes coincide with the eigenvectors of 



the instantaneous Hamiltonian at i = and t = tf. As 
discussed later in more detail the boundary conditions 
on the second derivatives, and consequently the condi- 
tions for <7o in Eq. (I13p are special, in the sense that we 
shall allow for discontinuities in the acceleration ijc at the 
edge times (in fact also elsewhere) . Physically this means 
that the trap is ideally allowed to be displaced suddenly 
a finite distance, inducing a sudden finite jump of the ac- 
celeration, whereas the velocity qc and the trajectory qc 
remain always continuous, qdt) can be interpolated by 
a simple polynomial ansatz that satisfies these boundary 
conditions. Once qc(t) is fixed, we get the trap trajec- 
tory qo{t) from Eq. ([3]). In principle there is no lower 
bound for tf [l3l |. However, there are always some lim- 
its in the laboratory related, for instance, to spatial or 
energy constraints. 



III. OPTIMAL CONTROL PROBLEM WITH 
CONSTRAINED RELATIVE DISPLACEMENT 

We begin with the equation of motion, Eq. ([3]), for the 
classical particle in the harmonic trap, and set, for com- 
pactness and to follow the usual conventions in optimal 
control theory, a new notation. 



xi = qc, X2 = 9c, u{t) 



qoi 



(14) 



where xi,X2 are the components of a "state vector" x, 
and the relative displacement between the trap and the 
center of mass u{t) is considered as the (scalar) control 
function. The physical motivation behind this control is 
that actual traps are not really harmonic, so the relative 
displacement should be kept bounded. Eq. ^ becomes 

Xl = X2, (15) 
X2 = -i^f)U. (16) 

The optimal control problem is to find \u{t)\ < 6 for 
some fixed bound 6, with u(0) = and u{tf) — such 
that the system starts at {a;i(G) — 0,a;2(0) = 0}, ends 
up at {xi{tf) — d,X2{tf) — 0}, and minimizes a cost 
function J. 

The boundary conditions for xi and X2 can be equiv- 
alently considered as those for qc and qc- The boundary 
conditions for u{t) are equivalent to those for q^ and, 
through Eq. equivalent to those for ijc, so there are 
totally six boundary conditions, as in Eqs. (|lll) and (|12p . 
A natural way to understand the boundary conditions on 
u{t) is to consider that u{t) = for t < and t > tf, 
so the center of mass and the trap center coincide before 
and after the transport. We will consider cost functions 
that are not affected by the isolated values u{0) and u{tf), 
for example minimizing the transport time or the energy, 
and solve the control problem in the interval {0,tf). In 
order to match the boundary conditions at the initial and 
final times, the optimal control obtained may be comple- 
mented by appropriate jumps at these points which do 
not affect the cost. We use Pontryagin's maximum prin- 
ciple, which provides necessary conditions for optimality 



3 



[38|. Generally, to minimize the cost function 



Jiu) 



g[yi{t),u]dt, 



(17) 



the maximum principle states that for the dynamical sys- 
tem X = f[x(t), u], the coordinates of the extremal vector 
x(t) and of the corresponding adjoint state p{t) formed 
by Lagrange multipliers, pi, p2, fulfill the Hamilton's 
equations for a control Hamiltonian He, 



x= ^ 
dp 

OH, 

P = 

OX 

where He is defined as 

He [p(<) , x(t) , u] = pog [x(t) ,u]+p'^ ■ f [x(0 , 



(18) 
(19) 

(20) 



The superscript "T" used here denotes the transpose of 
a vector, and po < can be chosen for convenience since 
it amounts to multiply the cost function by a constant. 
The (augmented) vector with components {po,Pi,P2) is 
nonzero and continuous. For almost all < t < tf 
the function Hc[p{t),x(t),u\ attains its maximum at 
u — u{t), and Hc[p{t),x{t),u(t)] = c, where c is con- 
stant. 



A. Time minimization 

We discuss now the time-minimization optimal control 
problem with a constrained relative displacement, that 
is, \u{t)\ = \qc — qo\ < ^, which means Ep < imwgJ^. To 
find the minimal time tf we define the cost function 



Jt 



dt = tf. 



The control Hamiltonian Hc[p{t),x{t),u] is 

Hc{pi,P2,Xi,X2,u) ^Pq +P1X2 -p2Wol 



With the control Hamiltonian, Eq. 
ing costate equations. 



(21) 



(22) 



Pi 

P2 



0, 



-Pl- 



191) gives the follow- 

(23) 
(24) 



They are solved easily aspi = ci andp2 = ~cit + C2 with 
constants ci and C2 . According to the Pontryagin's max- 
imum principle, the time-optimal control u{t) maximizes 
the control Hamiltonian in Eq. ((22|) . 

Since the control Hamiltonian is a linear function of the 
control function u{t), the optimal control that maximizes 
He is determined by the sign of p2, when u{t) is bounded, 
|u(i)| < When p2 0, the optimal control in the 
duration i/ is given by 



nit) = ^ T^' > 



(25) 
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FIG. 1. (Color online) (a) The control function (dashed blue 
line) for time-optimal problem, and (b) the corresponding tra- 
jectory (solid red line) for luq = 2n x 50 Hz, d = 1.6 mm, 
S = 0.16 mm and tf = 20 ms given by Eq. p8)) . 



If P2 = for some time interval, then pi = from Eq. 
(Ell), and Po = from Eq. since H^ = for the 

time optimal problem [38|, in contrast with the maxi- 
mum principle that requires (po,Pi,P2) 7^ 0. Thus p2 can 
be zero only at isolated points, the switching times. The 
solutions of the costate functions in Eqs. (1^ and ((M)) 
imply that the function of p2 depends linearly on time 
t, so that the sign of p2 cannot change more than once. 
Since the final point is {xi,X2) = {d,0),d > 0, the ap- 
propriate control sequence is of "bang-bang" (piecewise 
constant) type, 



L{t) = 



0, t<0 

-S, < t < ti 

S, ti <t <tf 

0, t>tf 



(26) 



with only one intermediate switching time at ti, as shown 
in Fig. [1] (a). The saturation of the control is typical of 
time minimization problems. 

Substituting u{t) into the classical Eq. ([2]), and using 
the boundary conditions in Eqs. (fTTj) and ([72]) . we find 
the optimal classical trajectory 



0, 



t < 
< t < ii 



u;lSty2, 

l^^^'-S d-Lol5{t-tff/2, ti<t<tf ' 
yd, t>tf 

and the corresponding trajectory for the harmonic trap 



4 



4 (^o) 




FIG. 2. Velocity profiles of harmonic trap and particles for 
time minimization, where the switching time ti = tf/2. 



qait) = 



0, t<0 

{l+uj^t^/2)S, 0<t<ti 

d- [ujl{t-tff/2 + l]6, ti<t<tf 

d, t>tf 



Fig. [T] (b) illustrates the time-optimal trajectory with 
one switching time. Solving the system of Eqs. dTSl) and 
(|16p. one can find the switching time ti and final time tf, 



(27) 
(28) 



by imposing continuity on xi and X2- For the "bang- 
bang" control the motion of the trap has discontinuities, 
while we impose continuity for the trajectory of the par- 
ticle. As illustrated by Fig. [51 the velocities of particle 
and trap become equal, 




qc = qo = 



0<t<ti 



ol5{t~tf), ti<t <tf ' 



(29) 



since u{t) is piecewise constant during the "bang-bang" 
control. The maximum velocity occurs at t = tf/2, 



(30) 



vq — uj^Stf /2 — (jJqV dS 



which is restricted by the imposed bound \u{t)\ < S. In 
addition, the instantaneous potential energy {V) is con- 
stant, and 



1 2r2 ^TTid^ 
Ep = -muj„S = -2-^. 



(31) 



If we loosen the bound by increasing S, the maximum 
velocity and the instantaneous potential energy increase, 
and the final time may be shortened. 



B. Displacement minimization 

In this subsection, we minimize the integral, or time- 
average of the relative displacement, which is equivalent 
to a minimal control-effort problem. To this end, the cost 
function can be defined as 



Jd 



\u{t)\dt 



\qc - qo\dt, 



/O .70 

and the control Hamiltonian is 



(32) 



(33) 



which leads to the same costate equations for pi{t) and 
P2{t) in Eqs. (|23p and We use for convenience the 

normalization po = —loq. Disregarding it-independent 
terms in He, the function of u(t) that we have to maxi- 
mize is 



P2U) = 



^ I -wg(l +P2)U, M > 
^ 'o(l-P2)u, M<0 



(34) 



According to Pontryagin's maximum principle, when u(t) 
is bounded, \u(t)\ < S, the control function is 



-S, P2 > 1 
u{t) ^ { 0, -1 <P2 < 1 
6, P2 < -1 



(35) 



which maximizes the control Hamiltonian in Eq. p3l) . 
Notice that whereas in the minimum-time problem dis- 
cussed above, the optimal control is "bang-bang" , the 
minimal-displacement control can be described as "bang- 
off-bang" , if we assume no singular intervals here. Ow- 
ing to the properties of costate equations, the "bang-off- 
bang" trajectory with two switching times ti and ti + t2 
can be described by, see Fig. [3] (a), 



0, t<0 

-S, <t<ti 

u{t) = ■( 0, ti < t < ti + t2 

S, ti+t2 <t <tf 

0, t>tf 



(36) 



Substituting the control function u{t) into the classical 
Eq. using the boundary conditions in Eqs. (fTTjl 

and (jl2p . and imposing the continuity of qc at the two 
switching times, the optimal trajectory for the center of 
mass, as shown in Fig. [3] (b), is finally given by 



qdt) = 



0, 

uji6ty2, 
vot-vl/{2ujlS), 
d-Lol5{t-tfY/2, 
d. 



t < 
0<t<ti 

h <t <ti+t2 

ti+t2 <t <tf 

t>tf 



which results in the following optimal trap trajectory, 

0, t<0 

{1 + uj^t^ /2)6, 0<t<ti 

qoit) ^ { Vat- vII{2uj15), ti < t < ti + t2 

d- [ujl{t-tf)^/2 + l]S, ti+t2<t<tf 

d, t>tf 
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FIG. 3. (Color online) (a) The control function (dashed 
blue line) for displacement-optimal problem, and (b) the 
corresponding trajectory (sohd red line) for i/ = 30 ms, 
S — 9d/2u!Qt'j, and other parameters are the same as Fig. 

El 



As in the "bang-bang" time-niinimization, the velocities 
of particle and trap are equal, see Fig. |4l 



vq, 



0<t<ti 

ti <t <ti+t2 



'UJ^S{t ~tf), ti+t2 <t <tf 



(37) 



where vq is the maximum velocity of trap motion in 
the trajectory, which will be determined later. With 
the boundary conditions for xi and X2 a-t t = ti and 
t ^ ti + 12, the switching times can be calculated as 



Vq 



d 

t2 = — 



Vo 



Vq L0^5' 

As a consequence, the final time is 



d vn 2 d 
2ti+t2^ — + ^> — a/t. 

uo WqO Wo V 



(38) 
(39) 

(40) 



Since the final time is fixed, there are three possible 
cases: (?) when tf > (2 / ljq) y^'d/S , the maximal velocity 
Vq can be solved from Eq. (HO)) as 



1 ± 



(41) 



where Vq should be ignored, because it leads to 2ti > tf. 
(ii) li tf — {2 / ujq) y/d/S , the maximum velocity is vq = 
u}Q\/dS, thus ti — tf/2 and t2 = 0. The trajectory in 



4 (^o) 




FIG. 4. Velocity profiles of harmonic trap and particles for 
displacement-minimization with two switching times ti and 

tl+t2. 



this case is reduced to that of the time-optimal control 
problem, (iii) When the time t/ is less than (2/ujQ)^/d/S, 
there is no real solution to vq and no solution to displace- 
ment minimization. 

Interestingly, the "bang-off-bang" trajectory obtained 
from displacement minimization may be related to the 
trajectory used for atomic transport in [1, 01 , where the 
shift velocity qq was increased linearly during a quarter of 
the spatial transported distance d/4, then kept constant 
for d/2, and during the last quarter finally ramped back 
to zero. To understand this in the context of optimal 
control theory for fixed tf and d, we note that for the 
choice 



9d 
2^f' 



(42) 



the maximum velocity is, according to Eq. (|4ip . 

M 



2tf' 



(43) 



and the switching times in Eqs. ((38)) and (p9| become 

(44) 



d _1 _ d 
2^' ^ ~ 3 ^ ~ 2^' 



The positions of the classical particle at the two switching 
times are 



9c(ii) = rf/4, qciti+t2) = M/A. 



(45) 



Due to the discontinuity at i = and t ~ tf, the motion 
of trap begins with grg(t = 0^") = 5 and ends up with 
qo{t = tj) = d — 6, so that 

qoit^) = d/A + 6, qo[iti+t2)+]=3d/A-S. (46) 

In other words, the protocol followed in 0,3 minimizes 
the averaged displacement by imposing the bound in Eq. 
(132) to the displacement. 



6 



Returning now to the general case, the time-averaged 
potential energy for the optimal trajectory is 



_ J^' Epdt muj^S'^ti 

tj-n — \ — \ : 



where ti = vq/lu'^S is given by 




As a result, 



-moj. 



2^2 



(47) 



(48) 



(49) 



For example, when 6 = 9d/2u!Qt^ and ti — tf/3 are cho- 
sen as discussed above, the time-averaged potential en- 
ergy is Ep = 27m(P /AuiQtp which is less than the (con- 
stant) potential energy Ep = SmcP/ujQt'j for the time- 
optimal control problem. 



C. Energy minimization 

The instantaneous potential energy {V{t)) is given in 
Eq. pop . To minimize the potential energy average for 
a given n and fixed transport time tf, the cost function 
can be defined as 



Je 



Epdt 



■f 1 



/o Jo 
and the control Hamiltonian is 
1 



—mcjjQU^dt, (50) 



He = -po-mujlu^ 



P1X2 



(51) 



which gives two costate equations, Eqs. and (1^^ . 

The solutions are pi = ci and p2 — —cit + C2, with 
constat ci and C2. For the normalization po = — 1/m the 
function of u{t) that we have to maximize is —v?/2—p2U. 

Here we start with the case of "unbounded control", 
i.e., without imposing any constraints on the displace- 
ment, and we shall show how this is related to the phys- 
ically interesting case where the control is bounded. To 
maximize — — P2U, the control function is found to 
be 



U{t) = -P2, 



and the classical Eq. g'c 
trajectory 



(52) 

-WqU, gives the optimal 



qc = -^ciult^ + ^02^0^^ + czt + C4. (53) 

Using the boundary conditions for qc and qc in Eqs. (jlip 
and (ini), we find ci = 12d/wg<^, C2 = 6d/wgiJ, cg = 



and C4 = 0. Clearly, Eq. (j53p docs not satisfy the bound- 
ary conditions for g'c in Eqs. (ITT|) and ((H]). To guarantee 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
t/tf 




FIG. 5. (Color online) (a) Control functions for the energy- 
optimal problem with bounded (solid red line) and unbounded 
(dashed blue line) controls, and (b) corresponding trajecto- 
ries in the cases of bounded (solid red line) and unbounded 
(dashed blue line) controls. The parameters are the same as 
in Fig. El 



u{t) = at i < and t > tf and match the boundary 
conditions, the control function u{t) has to be comple- 
mented by the appropriate jumps at these two edges. 
Consequently, the control function for unbounded con- 
trol, see Fig. [5] (a), is found to be 



0, 
0, 



t < 
0<t<tf 

t>tf 



(54) 



As shown in Fig. [5] (b), the optimal classical trajectory 
for unbounded control finally becomes 



t < 



|(3-2A), 0<t<t/ 



(55) 



t > t 



f 



where the trajectory qc in the interval (0,t/) is in agree- 
ment with the result obtained in using the Euler- 
Lagrange equation. In this case, the time-averaged min- 
imal potential energy is 



Ej, 



Jo 



Epdt 



6md 



(56) 



which gives a lower bound for the time averaged potential 
energy of any other trajectories satisfying all the bound- 
ary conditions, Ep > Gmd^/wgi^. Note that, in spite of 
not having preimposed a bound for the displacement, the 
optimal trajectory obeys \u{t)\ < Sq — Gd/utQt'j. For 
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FIG. 6. Velocity profiles of harmonic trap (dashed line) and 
center of mass (solid line) for energy-minimization with two 
switching times ti and ti + t2. 



the bounded control, i.e., when \u{t)\ < 6 is imposed, if 
5 > 5o the unbounded solution is the optimal one. (The 
value of 6o can be obtained in the bounded control case 
by requiring ii > 0, see Eq. below.) 

When the bound, \u{t)\ < S, is imposed, the control 
function is 



-5, P2> S 

u{t) = { -p2, -S <P2 < S 
S, P2 < -5 



(57) 



to achieve the maximum value of the control Hamiltonian 
He- As before, the linear p2 implies two switching times 
ti and ti + ^2- To make the control function continuous 
at ti and ti + ^2, it has the form shown in Fig. [S] (a). 




t < 
0<t<ti 

tf/2), ti<t<ti+t2 
ti+t2 <t <tf 
t>tf 



(58) 



where, because oitf — 2ti +t2 due to the symmetry, the 
two switching times ti and t2 are given by 



t2 = 26/ci, h = 



tf - 2(5/ci 



(59) 



Unlike the time-minimization and displacement- 
minimization problems, the control function here is not 
piecewise constant, so the velocities of the classical parti- 
cle and the trap are not equal during the second segment 
from ti to ti +t2, see Fig. [HI According to the control 
function in Eq. (|58p . imposing the boundary conditions 
for X2 at t = and t = tf, the velocity for the center of 
mass is 



LOQSt, 



-ulSit-tf), 



0<t<ti 
vo, h <t <ti+t2 ,(60) 

ti+t2 <t <tf 



and qo = qc — ii gives the velocity profile of the trap, 
46t, 0<t <ti 

-^ujQciit - ^-)'^ +vq- ci, ti < t < ti + 12 



qo = 



ti+t2<t <tf 



(61) 



where uq is the maximum velocity. With t2 = 25 /ci, and 
further imposing continuity of 2:2 at t = <i and t = ti+t2, 
we find 



vo_ 



_6_ 

2ci 



which finally leads to t/ = 2ti + 12, 

2t;o ^ 5 
*/ ^ ^ • 



(62) 



(63) 



Solving Eqs. \Q2\ and ([63|. the parameters ci and wq are 
given by 



Cl 



25 



1 



tf-2ti 4 



Vo = -LO^Sitf +2ti). 



(64) 



Thus, C2 = Stf /{tf — 2ti). So far, Ci, C2 and vq are 
all functions of ti. To determine ti we write down the 
optimal-energy classical trajectory from Eq. (j60p . 



r 0, 



qdt) 



2^0 
1, 

2 



t < 
< < < ti 

2; +Vot + C3, ti<t<ti+t2 



ti+t2 <t <tf 
t>tf 



By using the continuity of a;i at i = ti and t = ti + t2, 
C3 and ti can be solved as 



1 



C3 
ti 



{d - votf), 



(65) 

(66) 



where the other unphysical solution should be neglected. 
Once ti is fixed, cj (j = 1, 2, 3) are available, and vq is 
given by 



Vq = 



1 - 



%/3 



'1 - 



(67) 



which is less than the maximum velocity for the 
displacement-optimal trajectory. A trajectory with min- 
imal energy and bounded control is depicted in Fig. [5] 
(b). It is seen from Eqs. ([66| and (l67l) that for a real 
ti and Wo, tf > {2 / ujq) ^J'dJ5 should be satisfied. In the 
particular case tf — (2 / ljq) ^/dJS , the maximum veloc- 
ity is Vq — ujQ\/d6, thus ti — tf/2 and ^2 — 0. Like for 
displacement minimization, the trajectory in this case is 
reduced again to that of the time-optimal control prob- 
lem. Moreover, to make ti non- negative, tf should be 
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6000 
4000 



0.020 0.022 0.024 0.026 
tf {m$ 



0.028 0.030 



FIG. 7. Dependences of time-averaged energies Ep on tf with 
S — O.ld, for energy- minimization with bounded control (red 
soUd hne), unbounded control, Eq. (jSGjl (blue dashed line), 
and displacement-minimization (purple dotted line). Param- 
eters are d = 1.6 mm, a;o = Stt x 50 Hz, Eo = 100 x {hujo/2) 
and mass of Rubidium 87. 



less than {■s/&/ijjfj)yJ d/6. litf > {^/6/ujo)y/d/5, the opti- 
mal trajectory is the one in the unbounded-control case, 
as commented before. In other words, S > 6o = Qd/iJ^t^j. 
As a result, the segmented form in Eq. (j58p applies for 
the interval Ad/uiQt'j < S < Gd/uQt'j marked by vertical 
lines in Fig. [T] There is no solution for smaller times, 
whereas the solution becomes the one for unbounded con- 
trol for larger times. 

In this energy-optimal trajectory, the time-averaged 
potential energy Ep should be minimized. The cost func- 
tion in Eq. ([501 becomes 



Je = mwofS^ti -I- ^mujQS'^t2, 



and therefore 



lo' Epdt 



2c2 / 2ii 1 



which finally results in 



E.„ 



2 r2 
— mUJnO 



2V^ 



4d 



(68) 



(69) 



(70) 



In Fig. [71 we compare this to the (larger) average 
energy for the displacement-optimal problem, Eq. (I49L 
and the lower bound Eq. ([55)1 . and also demonstrate 
that the lower energy bound can be realized when t/ > 

(v^/l^o)v/^- 



IV. DISCUSSIONS AND CONCLUSIONS 

We have proposed optimal protocols for fast atomic 
transport in harmonic traps combining the invariant- 
based inverse engineering method and optimal control 



theory. Optimal trajectories with "bang-bang" and 
"bang-off-bang" forms are respectively obtained for time- 
minimization and displacement-minimization with con- 
strained displacement between the trap center and the 
center of mass of the particle density. The transient en- 
ergies for bounded and unbounded displacement are also 
minimized. 

In the time-optimal problem, the minimal time, Eq. 
([28|. corresponds to a fixed constraint 5. Consistently 
with this, no solutions are found for displacement and en- 
ergy minimization problems for transport times shorter 
than the minimal time, i.e. for tf < {2 / luq) y/dJS . To 
achieve fast and faithful transport in shorter times, an 
"energy price" must be paid by increasing i5 which, in real 
traps, will also produce errors because of anharmonici- 
ties. The relation between the minimal (time-averaged) 
energy and the transport time tf obtained here is not 
at all trivial, in particular they are not simply inversely 
proportional, see e.g. Eqs. ([5S1) or (ITD)) . as one might 
naively expect from the form of time-energy uncertainty 
relations. The scaling laws found are also peculiar of 
transport. For example the minimal energy in Eq. (j56p 
depends on t instead of the 1 7^ dependence applicable 
to engineered trap expansions [17| . 

In a previous work on invariants and transport [l3| , the 
energy bound for Ep was found using the Euler-Lagrange 
equation. Here we have shown how to realize this bound 
by allowing the discontinuous acceleration of the trap at 
t = Q and t = tf in the unbounded control optimization. 
In principle these and other discontinuities found could 
be avoided by imposing appropriate bounds and using a 
powerful pseudospectral numerical optimization method 
j33, 39, ,40] to address the corresponding more complex 
optimal control problem. 

Anharmonicity could be dealt with in a completely dif- 
ferent way using the protocols for anharmonic transport 
described in [l3|, which require a compensation of iner- 
tial forces in the frame of the trap. This may be feasible 
or not depending on the accelerations imparted and the 
corresponding optimization will be considered elsewhere. 

Last but not least, the present results may be extended 
to Bose-Einstein condensates following [ij]. Tonks- 
Girardeau gases could also be treated with a simple gen- 
eralization (16i] . 
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